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1 基于 16S rRNA 测序 分 析 奶 牛 体外 发 酵 中 瘤胃 细菌 多 样 性 的 变化 ! 
2 eee RH MOT 张 林 kak BIKAR” 
3 CPA IK SRS be, FARE 210095) 


4 W 要 : 本 试验 则 在 通过 16S rRNA 测序 技术 ， 分 析 奶 牛 体外 发 酵 瘤 胃 微 生物 区 系 随时 间 推 
立体 外 发 酵 分 析 方 法 提供 支持 。 试 验 使 用 体 
液 ， 与 人 工 唾液 混合 后 , 注入 盛 有 发 酵 底 


5 ” 移 的 变化 规律 , 为 探究 瘤胃 微生物 繁殖 规律 及 完 
6 ”外 产 气 法 进行 发 酵 。 采集 3 头 健康 奶牛 的 新 鲜 瘤 


TE 


7 WOES AER. 39 "C 水 浴 振 荡 培 养 ， 分 别 于 0、6、12、24、48h 采集 发 酵 液 ， 用 16S rRNA 


8 — 测序 技术 分 析 瘤 胃 细 菌 变化 规律 。 结 果 表 明 : D 所 有 的 序列 被 鉴定 为 14 个 门 和 59 个 属 。 


9 随 着 时 间 的 变化 ， 在 门 水 平 上 ， 纤 维 杆菌 门 〈(Fibrobacteres ) 、 拟 杆菌 门 (Bacteroidetes) 7I 


I. 10 — 胶 球 形 菌 门 〈Lentisphaerae) 的 相对 丰 度 减少 ， 纤 维 杆菌 门 、 拟 杆菌 门 在 48 h 时 显著 低 于 其 
= 1] ”他 时 间 点 (P<0.05)〉 ， 黏 胶 球 形 菌 门 在 24h 时 显著 低 于 其 他 时 间 点 CP<0.05) 。 而 螺旋 体 
e 12 [] CSpirochaetae) 和 变形 菌 门 (Proteobacteria) AN AA XTSE REINO, BREAST ITE 48 h 时 显著 高 于 
= 13. ”其 他 时 间 点 (P<0.05) ， 变 形 菌 门 在 12、24、48h 时 显著 高 于 其 他 时 间 点 CP<0.05) ; 随 
S 14 ”着 时 间 的 变化 ， 在 属 水 平 上 ， 丁 酸 弧 菌 属 CButyrivibrio) 和 瘤胃 球菌 属 (Ruminococcus) 的 相 


15 对 丰 度 增加 , 均 在 在 48 h 时 显著 高 于 其 他 时 间 点 (P<0.05); 而 解 政 珀 酸 菌 属 CSucciniclasticum ) 


16 ”和 普 雷 沃 氏 菌 属 (Prevotella) 的 相对 丰 度 则 减少 ， 解 政 珀 酸 菌 属 在 24、48 h 显著 低 于 其 他 时 


17. TR (P<0.05) ， 普 雷 沃 氏 菌 属 在 48 h 时 显著 低 于 其 他 时 间 点 (P<0.05〉 . 20 试验 获得 15 


18 ”792 个 操作 分 类 单位 (OTU) 并 进行 多 样 性 分 析 ， 稀 释 曲 线 达 到 平台 期 ，Simpson 指数 和 
19 ”Shannon 指数 均 表 明 菌 群 丰 度 具有 先 降低 后 升 高 的 趋势 。 综 上 所 述 ， 奶 牛 瘤胃 体外 发 酵 过 程 
20 F, 瘤胃 细菌 主要 的 门 以 及 属 的 相对 丰 度 以 及 多 样 性 均 随 时 间 推 移 发 生 改变 , HAE S 
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功能 差异 在 发 酵 过 程 的 变化 更 为 显著 。 


EE 


关键 词 ， 奶牛 ; 体外 发 酵 ， 微 生物 


的 功能 差异 与 底 物 以 及 人 了 
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因此 , 在 使 用 体外 发 酵 法 进行 研究 时 ， 应 充分 考虑 瘤 


[唾液 的 营养 物质 减少 引起 的 发 酵 环境 变化 的 关系 。 


中 图 分 类 号 : S823 


瘤 骨 是 一 个 复杂 的 生态 系统 ， 内 含 多 种 多 样 的 古 细菌 、 真 菌 及 原虫 等 ,它们 在 维持 瘤 


XA: 168 TRNA 


文献 标识 号 : 


文章 编号 : 


JI 


内 环境 稳定 和 影响 奶牛 健康 方面 起 关键 作用 。 这 些 瘤胃 微生物 具有 高 度 特异 性 , 它们 消化 代 
的 能 量 可 以 被 奶牛 利用 , 对 营养 物质 的 消化 起 着 关键 的 作用 , 还 可 能 会 提高 宿主 对 入 


谢 产生 


NAAR. 2 
质 主要 是 靠 纤 维 素 分 解 菌 ， 如 丝 杆 菌 (Fibrobacter) . 8 SEKE CRuminococcus) ~ T RINA 


C Butyrivibrio ) 等 B 


NP, X 


(Prevotella) Æ ARAH, "EM RI EA ERR 


菌 是 瘤 骨 中 数量 


最 多 的 一 种 微生物 ， 其 中 分 解 粗 饲料 中 的 纤维 物 


产 琥珀 酸 丝 杆 菌 具 有 很 强 的 分 解 纤维 素 能 力 ， 其 发 酵 产物 有 乙酸 和 琥珀 
! 菌 对 抗生素 也 有 很 强 的 耐 受 力 。 另外 ， 拟 杆菌 门 CBacteroidetes) 中 的 普 雷 沃 氏 菌 属 
白质 、 淀 粉 、 多 糖 ， 但 不 能 分 解 纤维 素 。 总 


之 ， 多 种 微生物 相互 作用 、 相 互 制约 ， 共 同 维持 瘤胃 内 环境 稳定 ， 对 于 促进 奶牛 消化 代谢 、 


维持 奶 
要 的 。 


产 气 法 是 应 用 较 广 、 操 作 较 简 
NE 


Fri a ASL TAS IHR AUR ARP, KE BAT CF DIS S 


R. dE PERSE, 


牛 健康 、 提 高 生产 性 能 具有 重要 意义 。 因 此 ， 对 瘤 骨 微生物 区 系 变化 的 研究 是 十 分 必 


发 酵 方 法 之 一 ， 常 用 来 评价 饲料 消 
究 底 物 精 粗 比 和 添加 剂 等 对 瘤 骨 微 


生物 区 系 的 影响 S73， 但 是 在 正常 饲 粮 作为 发 酵 底 物 时 ， 用 体外 发 酵 法 研究 瘤胃 微生物 区 系 


随时 间 变 化 的 报道 较 少 。 


下 的 科研 工作 具有 现实 意义 。 


16S rRNA 是 细菌 上 编码 rRNA 相对 应 的 DNA 序列， 存在 于 所 有 细菌 的 基因 组 中 ， 具 


有 高 度 


因此 ， 了 解 产 气 法 系统 带 来 的 瘤胃 微生物 区 系 变化 ， 完 善 现 有 条 件 


的 保守 性 和 特异 性 。 随 着 PCR 技术 的 发 展 ，16S rRNA 测序 已 经 成 为 分 析 微 生物 区 系 


Bl, Ae E 
变化 规 得 


1 材料 与 方法 


组 成 、 研 究 生 态 系统 多 样 性 的 重要 方法 ， 其 
在 使 用 16S rRNA 测序 技术 ， 分 析 奶 牛 体 
， 为 探究 瘤 骨 微生物 繁殖 规律 及 完善 体外 发 


准确 性 、 全 面 性 、 精 确 性 远 优 于 其 他 传统 的 技术 


外 发 酵 瘤 骨 微生物 区 系 随时 间 推 移 的 
酵 分 析 方 法 提供 支持 。 
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1.1 体外 发 酵 产 气 系统 建立 


发 酵 底 物 为 江苏 省 南京 西 岗 奶 牛 场 提供 的 全 混合 饲 粮 , 精 粗 比 为 50:50, 经 烘 干 粉碎 后 ， 
过 20 Hm, fH. ARERR CHORI WE 1 。 


dl 饲 粮 组 成 及 营养 水 平 〈 王 物质 基础 ) 


Table 1 Composition and nutrient levels of diets (DM basis) 
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% 
原料 Ingredients 含量 Content 营养 水 平 Nutrient levels? 含量 Content 
玉米 Corn 20.86 产 奶 净 能 NEL/(MJ/kg) 6.89 
XX Barley 1.89 粗 蛋 白质 CP 18.32 
豆粕 Soybean meal 12.49 粗 脂 肪 EE 4.33 
棉籽 Cottonseed 3.41 酸性 洗涤 纤维 ADF 22.81 
预 混 料 Primix! 2.70 中 性 洗涤 纤维 NDF 34.58 
玉米 酒糟 及 其 可 溶 物 Com 

6.24 $5 Ca 0.41 
DDGS 
小 麦 青贮 Wheat silage 20.75 we P 0.38 
Bis F= Alfalfa hay 12.13 
羊 草 Chinensis wildrye 4.04 
甜菜 粕 Beet meal 12.26 
湿 啤 酒糟 Brewers grain, wet 3.23 
合计 Total 100.00 


D 预 混 料 为 每 千克 饲 粮 提供 The premix provided the following per kg of the diet: VA 13.5 


KIU, VD33.78 KIU, VE 54 KIU, Cu 20 mg, Fe29 mg, Mn 26 mg. Zn 92.75 mg, 12.16 mg, 


Se 0.59 mg, Co 0.14 mg. 


2 产 奶 净 能 为 计算 值 ， 其 余 为 实测 值 。NEL was a calculated value, while others were 


measured values. 


发 酵 液 采用 新 鲜 瘤 胃液 与 人 工 唾液 混合 而 成 。 新 鲜 瘤 胃液 选 自 3 头 泌 乳 中 期 的 成 年 健康 


荷 斯 坦 奶牛 。 具 体 方法 为 : 在 晨 饲 2h 后 ， 通 过 瘤 骨 瘘管 接 取 足 量 的 瘤胃 液 ，4 层 纱布 过 滤 ， 


灌 入 充满 CO 并 经 过 预 热 达 39 "C 的 保温 瓶 ， 灌 满 后 把 瓶 口 立即 盖 紧 ， 送 回 实 验 室 ， 再 经 4 


层 纱 布 过 滤 ， 磁 力 搅拌 器 混 匀 ， 待 用 。 整 个 过 程 持 续 通 入 CO 维持 厌 氧 环境 。 人 工 唾液 参 


HE Theodorou 等 方法 配制 ， 人 工 唾液 组 成 见 表 2。 每 升 人 工 唾液 配置 过 程 为 : 在 558.9 mL 


蒸馏 水 中 加 入 0.1 mL 溶液 A. 200 mL 溶液 B、200 mL 溶液 C 及 1 mL RA D, ÑA CO» 
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饱和 后 放置 于 39 *C 并 恒温 水 浴 5~6 h， 加 入 40 mL 溶液 E， 混 匀 后 通 入 CO 至 饱和 并 加 热 


至 39*C， 继 续 通 CO, 0.$~1.0h， 使 人 工 唾液 变 为 无 色 ， 竺 用。 将 备 好 的 瘤胃 液 与 人 工 唾液 


在 厌 氧 条 件 下 充分 混合 (体积 比 1:9), pH 调 至 6.7，39 "CC 恒温 ， 作 为 发 酵 液 。 


表 2 人 工 唾液 组 成 


Table 2 Composition of artificial saliva 


项 目 溶液 体积 组 成 
Items Solution Volume/mL Composition 
A 微量 元 素 溶 液 100 CaCLb:2H;0 13.2 g, MnCb:4H?0 10.0 g, CoCl:6H2O 1.0 


g~ FeCl;6H;O 8.0 g 


B 缓冲 液 1 000 NH4HCO; 4.0 g, NaHCO; 35 g 

C 常量 元 素 溶液 1 000 NaoHPO, 5.7g、 KH2PO4 6.2 g、MgSO4.7H2O 0.6 g 
D JREF 1 0.1%《 质 量 体 积 比 ) 

E 还 原 剂 溶 液 100 1 mol/L NaOH 溶液 4mL、Na2S:9H2O 625.0 mg 


本 试验 使 用 160 mL RBA RA, BES A 1 g 发 酵 底 物 ， 持 续 通 入 
30s CO 后， 使 用 全 自动 分 装 仪 进行 快速 分 装 ， 每 瓶 100 mL 发 酵 液 ， 然 后 迅速 加 盖 异 丁 基 


塑胶 塞 ， 并 以 铝 盖 密 封 ， 摇 义 后 ， 置 于 39 "C 的 培养 箱 摇 床 中 ， 人 恒温 培养 48 h。 


1.2 ”试验 设计 


试验 采用 单 因素 试验 设计 。 待 培养 产 在 恒温 培养 0、6、12、24、48h 后 ， 分 别 取 出 放 


入 冰 水 浴 中 使 发 酵 停止 ， 并 分 别 吸取 2 mL 发 酵 液 置 于 冻 存 管 ， 密 封 后 液 氮 保存 。 本 试验 每 


个 时 间 点 设 3 个 重复 ， 共 15 个 样品 。 同 时 ， 另 取 3 个 培养 瓶 ， 只 加 入 人 工 唾液 进行 培养 ， 
用 于 设备 参数 校正 。 


1.3 ”累积 产 气量 的 测定 


培养 期 间 分 别 在 6. 12. 24. 48 h 测定 并 记录 体外 发 酵 产 气量 , 根据 朱伟 云 等 09 的 方法 ， 
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使 用 气压 转换 器 (OGER, UK) 5E 


时 测定 大 氧 瘤 


1.4 16S rRNA 测序 及 数据 分 析 


胃 微 生物 发 酵 产 气量 ， 根 据 各 产 气量 和 气压 
进行 校正 ， 除 去 空白 发 酵 瓶 产 气量 ， 计 算出 累积 产 气量 。 
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首先 对 样品 中 的 DNA 进行 提取 、PCR 扩 增 和 纯化 ， 然 后 构建 Miseq 文库 并 进行 测序 。 


接着 根据 index 和 barcode 序列 区 分 各 样本 数据 ， 且 对 数据 进行 处 理 ， 


优质 序列 长 度 分 布 。 利 用 QIME 软件 ， 使 用 RDP 分 类 器 L101 以 Silva、RDP、Greengene 细菌 


统计 有 效 数据 和 确定 


数据 库 呈 为 参 比 ， 对 序列 进行 系统 分 类 ， 分 类 基于 Bergey’s taxonomy， 共 分 为 6 级， 依次 


HF (kingdom) 、 门 (phylum) 、 纲 (class) 、 


BRUM EA 80%. Stir EAT FE i ZED) A 


gplots 软件 绘制 瘤胃 微生物 群落 热 图 "3。 最 后 ， 利 月 


位 (OTU) 69 进行 物种 多 样 性 分 析 。 


1.5 数据 统计 


采用 SPSS 17.0 软件 ANOVA 


示 ，P<0.05 为 差异 显著 。 


2 结 R 


2.1 测序 质量 与 合理 性 


过 程 进行 单 因子 方差 分 析 ， 试 验 数据 用 平均 值 + 标 准 差 表 


水 平 上 群落 组 成 ， 绘 制 系统 分 类 柱状 图 ， 并 采用 


H USEARCH 软件 031， 并 基于 操作 分 类 单 


H Corder) 、 科 (family) ~ J Cgenus) ， 


对 不 同样 本 DNA 进行 16S rRNA V3~V4 区 扩 增 ， 扩 增 的 条 带 大 小 正确 ， 浓 度 合适 且 条 


ii: 
dr 


带 单一 ， 说 明 纯 化 后 基因 组 DNA 


中 无 抑制 PCR 反应 的 杂质 ， 可 基于 16S rRNA 基因 序列 分 


析 微 生物 多 样 性 。 经 MiSeq 高 通 量 测序 平台 ， 共 得 到 748 172 条 raw. reads 序列 (每 个 样品 


平均 49 878 条 ) ; 经 去 除 序 列 两 端的 低 质量 区 域 、 其 中 的 接头 序列 、 含 有 模糊 碱 基 N 的 序 


列 程序 处 理 后 ， 对 序列 进行 组 装 拼接 ， 得 到 744 200 条 clean-reads 序列 (每 个 样品 平均 为 49 


613 条 ) ， 优 质 序列 比率 为 99.47%， 将 clean-reads 序列 去 除 引 物 错 配 过 多 的 序列 (>4bp) 、 


去 除 引物 序列 以 及 骨 合 体 后 ，15 个 样品 共 获 得 341 727 条 优化 序列 ， 每 个 样品 平均 22 782 


条 ， 其 


样品 稀释 曲线 是 采取 随机 抽检 


其 中 最 多 为 31 305 条 ， 最 少 为 14 195 条 ， 序 列 平均 长 度 456 bp. 


的 方法 获取 数据 ， 抽 到 的 序列 数 与 


所 代表 的 OTU 数目 
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构建 的 曲线 。 当 样品 稀释 曲线 趋 于 平缓 时 ， 说 明 测 序数 量 合理 ,更 多 的 测序 数据 对 于 发 现 新 


的 OTU 的 边际 效应 很 小 ， 反 之 ， 说 明 继 续 测 序 还 可 能 产生 较 多 新 的 OTU。 本 试验 15 个 样 


品 稀释 曲线 (图 1-A) 表明 并 非 所 有 样品 都 进入 了 平台 期 ( 即 随 着 抽取 数据 的 增加 仍 可 能 检 


测 到 新 的 OTU) ， 但 所 有 样品 的 履 盖 度 曲线 


(图 1-B) 都 已 达到 饱和 状态 ， 表 明 本 研究 样 


品 的 测序 量 及 测序 深度 合理 ， 更 深 的 测序 对 瘤 骨 液 中 检 出 新 的 OTU 贡献 很 小 。 


rarefaction curve 


> 


B coverage curve 


Number of OTUs 
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A: 样品 稀释 曲线 ; 


ed 


Coverge 
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B: hA. 


A: rarefaction curve; B: coverage curve. 


图 1 样品 稀释 曲线 和 覆盖 度 曲 线 


Fig. 1 Rarefaction curves and coverage - 


22 ”分 类 学 分 析 


从 图 2-A 可 以 看 出 ， 厚 壁 菌 门 (Firmicutes) 的 相对 丰 度 随 着 时 间 变 化 呈 先 上 升 后 下 降 的 


趋势 。 然 而 ， 拟 杆菌 门 、 纤 维 杆菌 门 CFibrobactere 的 相对 丰 度 总 体 上 呈现 下 降 的 趋势 。 螺 


旋 体 门 〈Spirochaetae ) 、 黏 胶 球 形 菌 门 CLentisphaerae) 和 浮 霉 菌 门 (Planctomycetes) 的 相对 丰 


度 呈 现 先 下 降 后 上 升 的 趋势 。 柔 膜 菌 门 (Tenericutes) 的 相对 丰 度 呈现 上 升 的 趋势 。 


由 表 3 可 知 , Act 48h 的 发 酵 ， 纤维 杆菌 
其 相对 丰 度 为 0.22%; 拟 杆菌 门 的 相对 丰 度 也 


门 的 相对 丰 度 显著 低 于 其 他 时 间 点 (P<0.05)， 


显著 低 于 其 他 时 间 点 CP«0.050 ， 其 相对 丰 度 


为 30.69%; 黏 胶 球 形 菌 门 在 24h 时 显著 低 于 其 他 时 间 点 CP«0.050 ， 其 相对 丰 度 为 3.41%。 
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浮 霉 菌 门 的 相对 丰 度 在 6、12、24P 时 显著 低 于 其 他 时 间 点 CP<0.05) , 24h 最 低 ， 其 相对 


丰 度 为 0.59%; JERE 


门 的 相对 丰 度 在 6h 时 显著 高 于 其 他 时 间 点 (P<0.05) ， 其 相对 丰 度 


为 36.69%; 变形 菌 门 的 相对 丰 度 在 12. 24. A8 h 时 显著 高 于 其 他 时 间 点 (P<0.05) , 24h 


最 高 , 其 相对 丰 度 为 18.23%; 螺旋 体 门 的 相对 丰 度 在 48 h 时 显著 高 于 其 他 时 间 点 (P<0.05)， 
其 相对 丰 度 为 7.05%; 柔 膜 菌 门 的 相对 丰 度 在 48 h 时 
对 丰 度 为 5.58%。 
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A: l'UE; B: 属 水 平 。 横 坐标 中 1 (1、2、3、4、5) 分 别 代表 0. 6. 12. 24, 48h. FBI 


3。 


A: phylum level; B: genus level. In the abscissa, 1 (1, 2, 3, 4, 5) mean 0, 6, 12, 24, 48 h, 


respectively. The same as Fig. 3. 


图 2 


在 体外 发 酵 条 件 下 瘤 骨 细菌 组 成 随时 间 变 化 分 类 柱状 图 


Fig.2 Classification histogram of rumen microbial composition with time under in vitro 


fermentation conditions 


A3 瘤 肯 微生物 在 体外 发 酵 条 件 下 细 戎 


组 成 及 其 相对 丰 度 ( 门 水 平 ) 


Table 3 Composition and relative abundance of rumen microbes under in vitro fermentation 
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conditions (phylum level) % 

项 目 时 间 Time/h P 
Items 0 6 12 24 48 SEM P-value 
拟 杆 菌 门 Bacteroidetes 41.10* 38.89^ 42.60* 37.31^ 30.69* 1.26 0.00 
厚 壁 菌 门 Firmicutes 22.15* 36.69% 21.43° 24.09° 29.10° 1.55 0.00 
黏 胶 球形 菌 门 Lentisphaerae 1426 8.93? 7.46° 3.41° 5.334 1.00 0.00 
Awil] Proteobacteria 3.12° 2.87° 16.75? 18.23? 16.32* 1.88 0.00 
纤维 杆菌 门 Fibrobacteres 6.93? 4.77° 2.214 2.90° 0.22° 0.62 0.00 
螺旋 体 门 Spirochaetae 2.97° 1.478 1.123 5.45? 7.05* 0.63 0.00 
FEA] Planctomycetes 1.40° 0.70° 0.66° 0.59° 1.148 0.10 0.00 
柔 膜 菌 门 Tenericutes 0.574 0.574 3.12* 4.69^ 5.58? 0.56 0.00 
TM7 0.44° 0.245 0.23° 0.20° 0.16° 0.03 0.01 
KEREI] Elusimicrobia 0.348 0.16° 0.16° 0.07* 0.15^ 0.02 0.00 
SRI 0.288 0.16° 0.165 0.115 0.08* 0.02 0.00 
WAIT Cyanobacteria 0.12? 0.12* 0.08? 0.06^ 0.07^ 0.01 0.02 
放 线 菌 门 Actinobacteria 0.07° 0.15* 0.095 0.07° 0.06° 0.01 0.02 
FF wl] Synergistetes 0.04 0.01 0.01 0.00 0.02 0.00 0.10 
绿 弯 菌 门 Chloroflexi 0.03* 0.01° 0.00° 0.01° 0.01° 0.00 0.00 

同行 数据 肩 标 不 同 字 母 表示 差异 显著 (P<0.05) ， 无 字母 或 相同 字母 表示 差异 不 显著 

(P>0.05) 。 下 表 同 。 
In the same row, values with different letter superscripts mean significant differenc 
(P<0.05) , while with no letter or the same letter superscripts mean no significant difference 
(P>0.05) . The same as below. 
从 属 水 平 细 菌 区 系 的 物种 相对 丰 度 分 析 ， 相 对 丰 度 在 0.01% 以 上 的 属 有 72 个 ， 不 能 归 
类 到 目前 已 知 菌 属 的 细菌 相对 丰 度 高 达 42.26%。， 从 图 2-B 可 以 看 出 ， 普 沃 氏 菌 属 
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CPevotella) 、 链 球菌 属 (Streptococcus) 和 琥珀 酸 弧 菌 属 CSuccinivibrio) 的 比例 随 着 时 间 


的 变化 呈现 先 上 升 后 下 降 的 趋势 ， 丁 酸 弧 菌 属 CButyrivibrio) 、 密 螺旋 体 菌 属 (Treponema) 


呈 先 下 降 后 上 升 的 趋势 ， 丝 杆菌 属 《Fibrobacter) 和 解 琥珀 酸 菌 属 (Succiniclasticum) 呈 下 


降 的 趋势 。 


由 表 4 可 知 ， 随 着 时 间 的 变化 ， 普 沃 氏 菌 属 在 12h 时 显著 高 于 其 他 时 间 点 〈P<0.05) ， 


其 相对 丰 度 为 31.19%; 链球 菌 属 的 相对 丰 度 在 6h 时 显著 高 于 其 他 时 间 点 CP«0.05) ， 其 相 


对 丰 度 为 23.03%; 琥珀 酸 弧 菌 属 的 相对 丰 度 在 12 h 时 显著 高 于 除 24h 外 的 其 他 时 间 点 


CP«0.05) ， 其 相对 丰 度 为 5.94%; RC9 gut group 的 相对 丰 度 在 24h 时 显著 低 于 其 他 时 间 


点 (P<0.05) ， 其 相对 丰 度 为 1.73%; 丁 酸 弧 菌 属 的 相对 丰 度 在 12 h 时 显著 低 于 其 他 时 间 点 


CP«0.05) ， 其 相对 丰 度 为 0.26%; 密 螺 旋 体 菌 属 的 相对 丰 度 在 12 h 时 显著 低 于 其 他 时 间 点 


| 


CP«0.05) ， 其 相对 丰 度 为 0.24%; 丝 杆菌 属 的 相对 丰 度 在 48 h 时 显著 低 于 其 他 时 间 点 


CP«0.05) ， 其 相对 丰 度 为 0.22%; 解 琥珀 酸 菌 属 的 相对 丰 度 在 24、48 h 时 显著 低 于 其 他 时 


间 点 CP«0.05D0 ， 其 相对 丰 度 分 别 为 2.22% 和 2.42%。 


表 4 瘤 骨 微 生物 在 体外 发 酵 条 件 下 细菌 组 成 及 其 相对 丰 度 ( 属 水 平 ) 


Table 4 Composition and relative abundance of rumen microbes under in vitro fermentation 


conditions (genera level) — 96 
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项 目 时 间 Time P fü 
Items Oh 6h 12h 24h 48h SEM P-value 
普 雷 沃 氏 菌 Prevotella 21.43° 28.26^ 31.198 28.79% 19.194 1.26 0.00 
解 琥珀 酸 菌 属 
431* 3.81° 3.92? 2.22° 2.42° 0.23 0.00 
Succiniclasticum 
RC9 gut group 3.16 2.938 2.36^ 1.73° 2.14° 0.15 0.00 
密 螺旋 体 菌 属 Treponema 1.15€ 0.54 0.24* 1.37^ 1.998 0.42 0.00 
丁 酸 弧 菌 属 Butyrivibrio 0.76? 0.48° 0.264 0.51% 0.978 0.07 0.00 
瘤胃 球菌 属 Ruminococcus 0.605 0.224 0.194 0.43* 0.852 0.07 0.00 
IHY IEEE Paraprevotella 0.528 0.41° 0.33% 0.224 0.31* 0.03 0.00 
厌 氧 弧 菌 属 Anaerovibrio sp. 0.52° 0.925 1.55* 1.77 1.122 0.16 0.05 
Candidatus saccharimonas 0.448 0.24? 0.23? 0.20? 0.16? 0.03 0.01 
琥珀 酸 弧 菌 属 Succinivibrio 0.36* 0.85* 5.943 5.28” 4.73° 0.64 0.00 
厌 氧 芽孢 杆菌 属 
0.11* 0.01^ 0.01^ 0.01° 0.02° 0.01 0.00 
Anaerosporobacter 
巨 球 型 菌 属 Megasphaera 0.01° 0.01° 0.06^ 0.32? 0.732 0.09 0.01 
链球 菌 属 Streptococcus 0.054 23.03* 7.25* 9.98? 8.29? 2.02 0.00 
甲烷 短 菌 
0.04° 0.02° 0.01° 0.02° 0.138 0.02 0.04 
Methanobrevibacter 
欧陆 森 氏 菌 属 Olsenella 0.03 0.01 0.01 0.04 0.04 0.00 0.21 
Ornatilinea 0.03* 0.01^ 0.00^ 0.01° 0.01° 0.00 0.00 
双 歧 杆菌 属 Bifidobacterium 0.01^ 0.128 0.078 0.02° 0.02° 0.01 0.00 
丝 杆菌 属 Fibrobacter 6.93? 4.77° 2.204 2.90° 0.22° 0.62 0.00 
Anaeroplasma 0.27 0.25 0.90 1.66 1.50 0.00 0.16 
165 
166 WE Cheat map) Ji FH Pure Ir E mr Eb Pte AE EO Yk HB Se AR, DL fi BA 
167 “地 聚合 大 量 数据 。 瘤 骨 内 18 种 门 水 平 的 细菌 进行 双向 聚 类 分 析 ， 结 果 显 示 (图 3-A) ， 可 
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从 横向 将 15 个 样品 分 成 2 个 大 簇 一 一 拟 杆菌 门 、 厚 壁 菌 门 ， 这 2 个 菌 门 聚 成 一 类 ， 其 余 16 
个 菌 门 聚 成 一 类 (又 可 到 成 2 7] RO ; 从 纵向 上 将 15 个 样品 分 成 2 DAG, 左边 的 第 1 大 
RASS 个 样品 (0，6 h〉， 右 边 的 第 2 KRAE 10 个 样品 (12、24、48h) 。 为 了 进一步 
了 解 瘤胃 微生物 细菌 相对 丰 度 随时 间 的 变化 规律 ,在 属 的 水 平 ， 随 着 时 间 的 变化 , 瘤胃 细菌 
组 成 存在 明显 差异 (图 3-B) ， 瘤 骨 内 59 种 菌 属 进行 双向 聚 类 分 析 结 果 显 示 ， 可 从 横向 将 
15 个 样品 分 成 3 个 大 簇 ， 普 雷 沃 氏 菌 属 聚 成 一 类 ， 链 球菌 属 聚 成 一 类 ， 其 他 菌 属 又 聚 成 一 
类 ;从 纵向 上 将 15 个 样品 分 成 2 DRE, AUR KRAE 3 个 样品 06h) ， 右 边 的 第 
2 大 簇 包含 12 个 样品 (又 可 聚 成 2 个 小 簇 ) 。 


A, Cale Key B Color Key 


Heatmap of phylum Heatmap of genus 


"e us Ila 


3 瘤 骨 微生物 在 体外 发 酵 条 件 下 细菌 的 双向 聚 类 分 析 


Fig.3 Two-way cluster analysis of rumen microbes under in vitro fermentation conditions 


2.3 ”测序 结果 及 a 多 样 性 分 析 


由 表 5 可 知 ， 随 时 间 变 化 ，Ace 指数 、Chao 指数 差异 不 显著 (P>0.05) ， 但 是 Shannon 
指数 与 Simpson 指数 差异 显著 (P<0.05)。 各 样本 稀释 曲线 均 达 到 平台 期 ， 说 明 测 序 量 及 测序 
深度 合理 ， 各 个 样品 的 覆盖 率 随 着 测序 的 深度 均 达 到 了 平台 期 ， 履 盖 率 都 在 97% 以 上 ， 说 
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184 。” 明 测序 数据 足够 大 ， 可 以 反映 样品 中 绝 大 多 数 微 生物 的 信息 。 


185 表 5 测序 结果 及 a 多 样 性 分 析 
186 Table 5 Sequencing result and a diversity analysis 
项 目 时 间 Time/h P (ü 
Items 0 6 12 24 48 SEM P-value 
操作 分 类 单元 数目 1 100.00 1007.33 1009.00 1 074.00 1 073.67 15.00 0.36 
OTU amount/ 个 
Chao 指数 Chao index 1 250.33 1 208.92 1 187.93 ] 241.88 1 222.65 10.66 0.80 
TT 
ta Ace 指数 Ace index 1 247.18 1 214.95 ] 200.99 ] 251.26 1235.71 9.89 0.80 
o " P : : 
N Shannon 指数 5.90? 4.83 5.29 5.21* 5.39 0.09 0.00 
e 
©) Shannon index 
qu Simpson 指数 0.01* 0.06? 0.06 0.02? 0.02? 0.00 0.00 
= Simpson index 
ON 覆盖 率 Coverage rate/% 98.00 98.00 98.00 98.00 98.00 0.00 0.66 
= 
x 187 2.4 pH 和 累积 产 气量 
e 
三 188 由 表 6 可 知 ， 随 时 间 变 化 ，pH 在 24 h 时 显著 低 于 其 他 时 间 点 (P<0.05) 。 随 时 间 变 化 ， 
O i89 ”累积 产 气量 逐渐 增加 ， 在 48h 时 显著 高 于 其 他 时 间 点 (P<0.05) 。 
190 表 6 pH 和 累积 产 气量 随时 间 的 变化 
191 Table6 Change of pH and accumulative gas production with time 
项 目 时 间 Time/h PE 
Items 0 6 12 24 48 SEM P-value 
pH 6.70* 6.68" 6.42» 6.26° 6.38° 0.05 0.00 
RRP AR 0.00° 80.534 172.40* 243.80^ 286.07" 23.59 0.00 


Accumulative GP 
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3 讨 论 


在 细菌 群落 中 


拟 杆菌 门 的 相对 比例 在 本 试验 中 是 最 高 


， 厚 壁 菌 门 、 
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拟 杆菌 门 和 变形 菌 门 是 瘤 骨 细菌 的 3 A OA a ES), e 


的 ， 这 与 之 前 给 美国 阿拉 斯 加 州 和 挪威 的 驼 鹿 、 梅 花 


鹿 、 水 牛 以 及 斯 瓦尔 巴特 群岛 的 驯鹿 饲 喂 冬 季 饲 粮 和 给 大 额 牛 饲 喂 树叶 和 稻草 的 混合 饲料 的 


DAZ RA BO), Qin 等 94 和 Tun 等 09] 也 说 


门 等 。 有 研究 表明 


瘤 骨 严格 厌 氧 的 环境 中 是 无 法 进行 生长 与 繁殖 。 拟 杆菌 
物 〈 的 主要 降解 者 ， 
起 主导 作用 。 因 此 ， 拟 杆菌 门 的 相对 


， 酸 杆菌 门 、 


FE 实 无 论 是 人 类 还 是 猎 科 胃 肠 道 菌 


菌 门 和 厚 壁 菌 门 均 占有 较 高 的 丰 度 。Tajima $20, Singh 等 RU 和 Peng 等 2 都 指出 荷 
和 水 牛 瘤胃 内 的 细菌 都 以 拟 杆菌 门 为 主 ， 其 余 菌 门 的 相对 丰 度 比较 少 ， 比 如 酸 杆菌 
蓝藻 门 的 细菌 更 适合 在 江海 湖泊 等 环境 中 生存 PI， 它 们 在 


门 是 瘤 骨 微生物 中 非 纤 维 性 碳水 化 合 


含有 与 降解 非 纤维 类 多 糖 有 关 的 基因 P9， 而 纤维 分 解 菌 在 发 酵 中 后 期 
FE 度 在 12 h 之 后 减少 。 奶 牛 瘤 骨 内 不 会 产生 纤维 降解 


酶 分 解 摄 入 的 纤维 物质 ， 而 是 依靠 体内 的 细菌 、 真 菌 等 微生物 对 其 进行 降解 消化 


而 本 试验 纤维 杆菌 


纤维 杆菌 门 消 化 和 分 解 纤维 物质 起 着 重要 的 作 月 


Ho Liu 等 29 的 研究 表明 ， 纤 维 杆菌 


杆菌 属 从 0 到 16 和 48h 的 相对 丰 度 显著 增加 ,在 降解 植物 纤维 性 物质 方面 发 挥 关键 的 作用 。 


群 中 , 拟 杆 
斯 坦 牛 


rl, ie 


。 其 中 
门 的 纤维 


门 在 12 到 48 h 降解 纤维 时 显著 减少 ， 这 与 发 酵 方法 不 同 有 关 。 本 试验 是 


采用 的 体外 发 酵 法 ， 随 着 时 间 的 变化 ， 上 共有 分 解 纤维 功能 的 细菌 之 间 会 相互 竞争 ， 纤 维 素 、 


半 纤 维 素 物质 会 不 断 减 少 , 这 将 会 导致 纤维 杆菌 


门 相 对 丰 度 显著 减少 。 而 Liu 等 29 采 用 尼龙 


袋 法 ,在 每 个 时 间 点 都 会 供应 相同 的 饲料 ， 这 样 可 能 会 组 解 纤维 分 解 菌 之 间 的 竞争 。 同 时 在 
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Abstract: This study was aimed to assess the dynamics of microbial bacterial community in the 
rumen of dairy cows with time in vitro fermentation analyzed by 16S rRNA sequencing 
technology, which provided a support for exploring the rumen microbial propagation and 
improving the in vitro fermentation method. The fermentation was carried out using an in vitro gas 
production method. The rumen liquid from 3 healthy Holstein dairy cows were collected, mixed 
with artificial saliva, and then injected into fermentation flasks containing fermentation substrates, 
and incubated at 39 °C. Samples were collected at 0, 6, 12, 24, 48 h. The dynamics of microbial 
bacterial community in the rumen was analyzed by 16S rRNA sequencing technology. The results 
showed as follows: 1) a total of 14 phyla and 59 genera were identified from all samples. With the 
change of time, at the phylum level, the relative abundance of Fibrobacteres, Bacteroidetes and 
Lentisphaerae declined, thatof Fibrobacteres and Bacteroidetes at 48 h was significantly lower 
than that at other time points (P<0.05), and that of Lentisphaerae at 24 h was significantly lower 
than that at other time points (P<0.05). The relative abundance of Spirochaetae and Proteobacteria 
increased, that of Spirochaetae at 48 h was significantly higher than that at other time points 
(P«0.05), and that of Proteobacteria at 12, 24 ,48 h was significantly higher than that at other time 
points (P«0.05). At the genus level, the relative abundance of Butyrivibrio and Ruminococcus 
increased, that of both them at 48 h was significantly higher than that at other time points (P<0.05); 
the relative abundance of Succiniclasticum and Prevotella decreased, that of Succiniclasticum at 
24, 48 h was significantly lower than that at other time points (P«0.05), and that of Prevotella at 
48 h was significantly lower than that at other time points (P«0.05). 2) The experiments were 
obtained 15 792 OTUS, diversity analysis was carried out, and the dilution curve reached the 
plateau stage. Simpson index and Shannon index showed that the abundance of bacteria first 
decreased and then increased. The results showed that the relative abundance and diversity of the 
main phylum and genus of rumen bacteria changed with time and the changes in fermentation 


process were more significant due to differences in rumen bacterial function. Therefore, when 
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using the in vitro fermentation method for research, the relationship between the functional 
difference of rumen bacteria and the fermentation environment caused by the nutrients of artificial 


saliva and the reduction of substrate should be fully considered. 
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"Corresponding author, professor, E-mail: zyhan6708(a)njau.edu.cn (责任 编辑 陈 XD 


